Sum Rule Anomaly from Suppression of Inelastic Scattering in the Superconducting 

State 
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In the conventional BCS description of a superconductor the kinetic energy increases in the 
superconducting state. We describe the observed decrease in kinetic energy by adopting a simple 
model of electrons whose elastic scattering rate undergoes a sharp decrease as the temperature 
is lowered below T c . This phenomenology has been suggested by other experiments, particularly 
microwave conductivity. We find that such a decrease accounts for the observed increase; a study 
of these different phenomena over a wide range of high T c materials would confirm this correlation. 

PACS numbers: 



I. INTRODUCTION 

Much of the focus in the high temperature cuprates has 
been on the normal state. In fact, although we continue 
to probe the superconducting state with finer probes and 
better samples, the general conclusion seems to be that 
the state is BCS-like with well-defined quasiparticles. In- 
dications of this statement are from photoemissioni and 
thermal conductivity 2 studies on Bi2212 and YBCO, re- 
spectively. These and many other studies are converg- 
ing on the general notion that, while the superconduct- 
ing state has many conventional properties (albeit with 
an order parameter with d-wave symmetry), the normal 
state is in many ways anomalous. 

However, optical studies over the last half dozen years 
have found an anomalous behaviour in the optical sum 
rule^, that appears in the superconducting stat o 4 : 5 : 6 . 
This sum rule relates the weight of the entire oscillator 
strength to the bare plasma frequency; a related sum 
rule pertains to a single band; then optical processes 
involving transitions of excitations that belong only to 
this band must be well-separated (even in principle) from 
those that involve other bands. Alternatively, apparent 
violations may occur because higher energy processes are 
"left off the accounting sheet" , so to speak. 

The observations, whose interpretation remains under 
debatejiSiS can be summarized as follows. First, recall 
the "single band" sum rule—l^H 
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where efc is the tight-binding dispersion and is the sin- 
gle spin momentum distribution function. We define the 
quantity in braces as — < T xx >. Note that in general 
the kinetic energy is given by < K >= ■£> ^ fc ekUk, and 
only in the case of nearest neighbour hopping in two di- 
mensions is < T xx >=< K >. Nonetheless we will use 
this case to build intuition concerning the behaviour of 
the sum rule. As explored in Ref. Il-'l this characteri- 
zation is reasonably accurate even when beyond nearest 
neighbour hopping is included. 



The simplest situation in which one can understand 
Eq. n is the case of non-interacting electrons; then, the 
momentum distribution is simply given by the Fermi- 
Dirac function. At zero temperature this is a step func- 
tion at the Fermi energy, and all the states with the low- 
est energy are occupied, while those with higher kinetic 
energy are empty. This, of course minimizes the total ki- 
netic energy, and therefore maximizes the sum rule. As 
we raise the temperature, states with higher kinetic en- 
ergy become partially occupied, while those with lower 
kinetic energy become partially empty. This behaviour 
is illustrated in Fig. (la). Therefore the RHS of Eq. 
decreases as temperature increases. The sum rule (as de- 
termined by the RHS of Eq^l is plotted in Fig. 2 (upper- 
most red, solid curve) as a function of T 2 . The result is 
mostly linear in T 2 , as one can discern from the Sommer- 
feld expansion 14 . However, note that this will depend in 
general on the details of the band structure, filling, etc.; 
careful inspection of Fig. (2) shows deviations from T 2 
behaviour. This is partly caused here by the two dimen- 
sional van Hove singularity. However, other band fillings, 
as we have verified, also cause significant deviation from 
T 2 behaviour. 

What happens when the system goes superconducting? 
The momentum distribution function is shown in Fig. 
(lb), for a gap of varying size in the quasiparticle spec- 
trum. Keep in mind that for an order parameter with d- 
wave symmetry, the momentum distribution is no longer 
a function of the band structure energy, ej; alone. For 
example, for a BCS order parameter with simple nearest 
neighbour form, = y( cos k x —cos k y ), as k varies from 
the (0, 0) to (it, 0), the magnitude of the order parameter 
changes from zero to A, On the other hand, as k varies 
along the diagonal (from the bottom of the band to the 
top), the order parameter is zero (and constant). In any 
event, Fig. (lb) conveys the well-known fact that even 
at zero temperature, BCS-like superconductivity raises 
the kinetic energy of the electrons. The consequence of 
this for the electron kinetic energy is shown in Fig. (2), 
for both an s-wave (green, dashed curve), and a d-wave 
(blue, dotted curve) order parameter. As our discussion 
indicated, the steady increase as the temperature is low- 



1.0 



|= 0.5 



0.0 









A = 


\\ 


A = 0.1 


\ \ 
\\ 
\\ 

v 


A = 0.2 


T = 


\ 

i\ 
\\ 
\ \ 




\ \ 
\ \ 




N ^ - 




_z 



c 0.5 




c 0.5 




FIG. 1: A summary of various reasons for why the momentum 
distribution is broadened. In (a) temperature is responsible 
for the broadening, while in (b) the occurrence of supercon- 
ductivity results in broadening, even at T = 0. In (c) and (d) 
inelastic scattering (as modelled by coupling to an Einstein 
phonon within the Migdal approximation) and elastic scat- 
tering (as modelled by impurity scattering in the Born limit) 
lead to the broadened distributions as shown. 
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FIG. 2: The optical sum rule (RHS of Eq. Q) as a function of 
temperature for the clean limit (upper curves) and with elastic 
scattering (lower curves). Note that elastic scattering mainly 
just shifts the curves. An s-wave order parameter (green, 
dashed curve) leads to a more significant decrease in the sum 
rule than a d-wave order parameter (blue, dotted curve). We 
used a tight-binding model with nearest neighbour hopping 
only, as described in the text. 



ered is arrested at T c , and the magnitude of the kinetic 
energy in general decreases below T c . As expected, the 
effect is stronger for s-wave than for d-wave symmetry. 

Briefly, then, the experimental situation^ is that, par- 
ticularly for underdoped and optimally doped samples, 
the conductivity sum rule (i.e. the left-hand side (LHS) 
of Eq. ^| indicates an increase in the magnitude of the 
kinetic energy (RHS) in the superconducting state. For 
the purposes of this work, we assume this is the " experi- 
mental fact" and defer doubts and debates to Refs. 
(a more recent work is Ref. Ull) . 

Interestingly, this very scenario was predicted^ well 
before experiments with the required accuracy were per- 
formed. The basis for this prediction was that the elec- 
tron (or hole) kinetic energy can, in principle, change 
upon entering the superconducting state, so that the 
"single band" sum rule, Eq. ^ would not remain con- 
stant. In particular spectral area from relatively very 
high energy would change at the superconducting tran- 
sition. Since this spectral weight is not accounted for 
in the single band sum rule of Eq. ^ then there would 
be an apparent violation. This was demonstrated to oc- 
cur for the hole mechanismiii& of superconductivity. In 
that model the effective mass of the holes decreases in 
the superconducting state; by pairing the holes are 'lib- 
erated' compared to their mobility in the normal state. 
Partly this physics manifests itself in an effective pairing 
interaction which has the peculiar form obtained in Ref. 
ITsl In this sense these experimental results represent a 
" smoking gun" that point towards the hole mechanism of 
superconductivity. Many details are worked out in Ref. 
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However, other explanations have been put forth after 
the experiments were performed. Many have focused on 
the impact of superconductivity on the bosonic spectral 
function, i.e. on the boson that is thought to mediate the 
pairing. For example, in Ref. [12| a frequency-dependent 
scattering rate was assumed, that underwent a significant 
reduction in the superconducting state. A similar mod- 
elling was implemented in Ref. [l9| and Ref. |2fJ. In each 
case, a change was put in by hand; in the first case they 
were largely motivated by the Angular Resolved PhotoE- 
mission Spectroscopy (ARPES) results^, while in the 
latter two the appearance of a sharp neutron resonance 
modeaa provided the motivation for an abrupt change. 
In fact, in Ref. [2(3 this latter scenario was explored in 
more detail with models used for the spin fluctuations 
that have described a number of other experiments. 

The purpose of this paper is to see to what extent 
changes to a frequency-independent scattering rate can 
reproduce the observed sum-rule variation. The motiva- 
tion is two-fold. First, it is always helpful to simplify the 
description of the phenomenon as much as possible, so 
that presumably only the essence remains. Secondly, this 
follows the microwave analysis of Ref. I23L and so would 
bring together observations from different experiments. 
In that work (see also Ref. l24f) a two-fluid interpreta- 
tion of their results led them to a frequency-independent 
scattering rate whose value collapses in the superconduct- 
ing state, following a T 4 power law at low temperatures. 
This properly accounted for the broad low temperature 
peak in the real part of the microwave conductivity. In- 
deed, the original microwave experiments 25 interpreted 
the presence of this peak in the same way; both regarded 
the origin of this lifetime decrease as being due to a " gap- 
ping" of some collective mode responsible for the normal 
state scattering, but no longer functional in the super- 
conducting state. This is somewhat consistent with the 
scenario adopted by Norman and PepiniS and Knigavko 
et ali^. Here, we wish to simplify the phenomenology 
that relates all of these experimental findings: a sim- 
ple number characterizes the single electron scattering in 
the normal state (here this really is a "characterization" , 
as the normal state is very likely not even a Fermi Liq- 
uid). Thus linewidths are in general broad, as seen by 
ARPES, for example, and collective modes of these elec- 
trons will also be broad. A sharpening of the electron 
spectral function in the superconducting state leads to 
quasiparticle-like peaks in the ARPES, a collective mode 
that also sharpens, the prominent peak in the microwave 
conductivity, and an apparent sum-rule anomaly below 
T c . We will focus on a description of the latter phe- 
nomenon. We start with a discussion of the momentum 
distribution function in the next section, where, in many 
cases an analytical form can be obtained. We then ex- 
plore the impact this has on the right-hand side of the 
sum rule equation, Eq. Q), as a function of temper- 
ature, order parameter symmetry, band structure, and 
clastic scattering rate. Using the phenomenological tem- 
perature dependence of the elastic scattering rate derived 



from the microwave measurements we then illustrate the 
"anomalous sum rule behaviour", as seen in experiment. 
We close with a summary and implications for experi- 
ments. 



II. MOMENTUM DISTRIBUTION FUNCTIONS 

Returning to Fig. 1, we outline our calculations and 
explain the other parts. In general, the momentum dis- 
tribution function is given by 

Uk = J2 G(k,iu m )e l ^ 0+ = - + - J2 ReG(fc,z^ m ), 

" rn— — 00 m— — 00 

where G(fc, iu) m ) is the (in principle) exact single electron 
Green function with momentum k and frequency ito m . 
The iuj m are the Fermionic Matsubara frequencies, iui m = 
TrT(2m — 1), where m is an integer and T is the temper- 
ature (ks = 1)- The inverse temperature, (3 = l/T. The 
Green function G is generally given by the Dyson equa- 
tion, 

G{k,z) = l/[z-e k -E(k,z)], (3) 

where z is anywhere in the complex plane, is the non- 
interacting dispersion, and z) is the electron self en- 
ergy. For the superconducting state, this equation holds 
when G is understood to be the Gn element of the 2x2 
Nambu matri5*S&. 

For the RHS of Eq. the calculation of the self energy 
for a given model and approximation suffices, and these 
equations are all that is required. For the LHS, the two- 
particle response is required, and is much more difficult. 
Although we have computed the LHS for various models 
in the past, this calculation is on less sure footing, so we 
will show calculations of the RHS only. First, for the non- 
interacting electron gas, Go(k,z) = l/[z — ejJ, and the 
result from Eq. [21 is the Fermi function. This is what is 
plotted in Fig. la for the various temperatures indicated, 
all in units where the bandwidth is 2, and the chemical 
potential (hidden in Ck) is in the middle of the band. 

In Fig. (lb) we require the superconducting state. In 
this case the relevant Nambu matrix element is 

G (k iu> ) - iu m Z(k,iu m ) + e fc 

(iuj m ) 2 Z 2 (k,iuj m ) ~ 4> 2 {k, iuj m ) - e\ ' 

(4) 

where the self energy is given by iu} m ) = iu> m \l — 
Z(k, iu> m )] , and an off-diagonal self energy, <p 2 (k, iu> m ) = 
Z{k,iu)m)l±.(k,iuj m ), is also required. In BCS-like the- 
ories it is the A(/c, iuo m ) that reduces to the familiar 
gap parameter, A. For the purposes of Fig. (lb) we 
have assumed that something has given rise to super- 
conductivity, with a BCS order parameter, A^., whose 
values (constant in k for this figure) are indicated. Then 
Z(k, iuj m ) = 1 and A(k, iu) m ) -> A fc . Using Eq. g| Eq. H 
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yields the familiar 26 



E k 



(5) 



where = \/e| + A|. A finite and constant order pa- 
rameter clearly gives rise to smearing at what was the 
Fermi surface. Note that this is in the clean limit, and 
at zero temperature. Nonetheless, the BCS ground state 
is usually considered a Fermi Liquid, in the sense that 
well-defined quasiparticles (with energy Ek) exist, even 
though there is no longer a discontinuity at the chemi- 
cal potential. It is clear that the onset of superconduc- 
tivity causes occupation of higher kinetic energy states, 
just like raising the temperature does. If the onset of 
superconductivity more than compensates the lowering 
of temperature, than the trend which the kinetic energy 
was following with temperature will reverse itself. Fig. 
2 clearly indicates that this is the case, though less so 
for an order parameter with d-wave symmetry than one 
with s-wave symmetry. 

In what other ways c an higher kinetic energy states be 
occupied? Refs. Il2lll9l used a coupling to a phonon-likc 
boson. Both applied the so-called Migdal approximation. 
Using a single mode at Einstein frequency u>e, coupled 
to the electrons with coupling constant A, the self energy 
is given by 



Y,i n (iuj m ) = — iAdj^tan 



(6) 



Then the momentum distribution function for inelastic 
scattering becomes 



1 



_ v 

/?tr e2 fc + K+ A ^ tan_1 (^)) 



(7) 



In the zero temperature limit, u) m — > u>, a continuous 
variable and the Matsubara sum becomes a continuous 
integral. Then it is easy to see that a discontinuity re- 
mains at the Fermi level, whose size is 1/(1 + A). Calcula- 
tions are shown in Fig. (lc). Clearly increasing the cou- 
pling tends to push electrons into higher kinetic energy 
states. If a lowering of the electron-boson coupling were 
to accompany the onset of superconductivity as the tem- 
perature is lowered, the savings in kinetic energy could 
more than offset the squandering of kinetic energy seen 
in Fig. lb due to the superconductivity itself. This has 
been amply discussed in Refs. ll2TH^ . and will not be 
further addressed here. 

The focus of this work is elastic scattering, for which 
we will adopt the simple Born approximation. This will 



make everything we compute look formally like impu- 
rity scattering (in the Born limit). However, as Hosseini 
et al^ pointed out, the amount of impurities in their 
samples is small, and this issue has been tested under 
controlled impurity addition, so we are really modelling 
some inelastic process which, as discussed in the intro- 
duction, seems to disappear when the material goes su- 
perconducting. Elastic scattering is like static inelastic 
scattering, i.e. in Eq. (jSJ uje — > while Xlue — > l/( 7rr ), 
with 1/t = r some characteristic scattering rate. Then 
T, e \(k,iuj m ) = — i^sgn(uj m ), as expected for the Born ap- 
proximation. Substitution into Eq. |2 then yields, for 
elastic scattering, 



n k,el 



1. 



-Im ip 



r 

2 



(8) 



where ip(z) is the diagamma function, and q = q/(2irT). 
For T = one recovers the simple result 

n k , e i = - cot" 1 (^) (T = 0). (9) 

7T 1 

This quantity is plotted in Fig. (Id). Note that at finite 
temperature (shown below) more impurity scattering re- 
inforces what a higher temperature already does, i.e. it 
increases the kinetic energy. Thus, it should be clear that 
a collapse of the scattering rate will decrease the kinetic 
energy, and offset the trend evident in Fig. (lb). 

Finally, we wish to examine the distribution function 
in the superconducting state. The quantities 



Z(k,iw m ) = 1 + 



1 



4>(k,iuj m ) = A k + 



1 A fc 



(10) 



are entered into Eq. J3J, which is then inserted into Eq. 
©. The result is 



1T-k,els = n k,els + Snk,els 



where 



1 £k 2 T hf 1 ^^-* 

n k,els = 7T - "F V o + 77 + lE k 

Z rj k 7T \ZZ 



(11) 



(12) 



just like Eq. (jHJ) except that E k replaces e k at relevant 
places, and the additional piece, 



J 



5n kMs = -- ^ 



+ (v^TAf + l/(2r)) 2 4 + K + V(2t)) 2 + A 2 fe 



(13) 
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has to be evaluated numerically. Nonetheless, in the most 
extreme case we looked at it represented less than 0.1% 
of the total result and can thus be safely ignored (we 
have included it in all our displayed results). Eq. H12|) 
is a nice compact form that is also particularly simple at 
zero temperature: 



1 A e k 2 
n k,ei s = r 1 - -7; — tan 

2 V ti k 7T 



(14) 



This form makes it clear that elastic scattering and the 
onset of superconductivity both tend to broaden the dis- 
tribution function, as indicated by Fig. 1. 



III. THE RHS 

Results for the RHS of Eq. (JTJ are most easily deter- 
mined by adopting a model density of states which is flat 
and finite; one can then perform the integral straight- 
forwardly numerically (and in some cases analytically). 
However, this can't be done when an order parameter 
with d-wave symmetry is present. Hence, we have used a 
tight-binding band structure, and simply performed the 
momentum sums (in two dimensions) in Eq. (T|. 

We begin with nearest-neighbour hopping only. Then 
results for the clean limit were already displayed in Fig. 
2. The additional curves are shown for an elastic scatter- 
ing rate 1/r = 5 meV. The dashed, green (dotted, blue) 
curves are for s-wave (d-wave) order parameters, respec- 
tively. Note that elastic scattering does almost nothing to 
the results except for a constant shift downwards, consis- 
tent with Fig. (Id). Even in the superconducting state, 
the presence of elastic scattering simply gives rise to a 
constant shift downwards. 

To what extent does the "standard BCS" picture 
shown in Fig. 2 depend on the use of a tight binding 
model with nearest neighbour hopping only? In particu- 
lar, the hole- (and electron-) doped cuprates have band 
structures that are often modelled with a next-nearest 
neighbour contribution to the hopping. To explore this 
issue in a little more depth, we define the relative dif- 
ference between the superconducting and normal state 
value of the sum rule as: 

< T _ (-<T XK>1 )- (-<T xx>jv ) 

5T t = , (15) 

(— < ixx >N) 

where i = S, D refers to the superconducting state with 
s- or d-wave symmetry, respectively, and the subscript 
N refers to the normal state, all at a given tempera- 
ture. We use a similar definition for the difference in 
the magnitude of the kinetic energy, 5{—Ki). Inspection 
of Fig. 2 (where 8Ki = 8Tj) indicates that this quan- 
tity is negative, with a value of about 1% for s-wave and 
about half that value for d-wave. The important point is 
that it is negative. To see the sorts of changes expected 
from a next-nearest neighbour amplitude and/or change 
of electron density, we show in Fig. 3 a plot of STi vs. 



t' (curves), at zero temperature, for a band structure, 
e.k = — 2£(cos k x a + cos k y a) + At' cos k x a cos k y a. We use 
t = 62.5 meV, and allow r = t' /t to range from —0.3 to 
0.3, for two electron densities, n = 1, and n = 0.85. We 
have taken T c = 80 K, and used BCS values for the or- 
der parameter as a function of temperature; for s-wave, 
2A /k B T c = 3.5, while for d-wave 2A /fc B T c = 4.2. For 
the temperature dependence we have solved BCS equa- 
tions in either case, with a small attractive on-site poten- 
tial for s-wave symmetry, and a small nearest neighbour 
attractive potential for d-wave symmetry, on a square 
lattice. Actually, we found that the temperature depen- 
dence is nearly the same i n both cases, a nd very well 
described by A(T)/A = y/l - (i 4 + i 3 )/2, over the en- 
tire temperature range, where t = T/T c . 

Fig. 3 shows that there are minor variations with fill- 
ing and band structure. We have also shown the rela- 
tive change in kinetic energy (symbols) to illustrate that 
these are qualitatively the same as the relative change 
in the sum rule, thus confirming van der Marel et al.'s 
observation^. Note that for electron doping one acquires 
the same results as for hole doping for negative t'. An- 
other trend to note is that the relative value of the effect 
due to superconductivity goes down as the overall scale 
of the electronic band structure increases with respect to 
the superconductivity energy scale (80 K ss 7 meV). 



IV. REVERSAL OF KINETIC ENERGY 
CHANGE AT T c 

Inspection of Fig. 2 (compare two red, solid curves) 
makes it clear that a sudden drop in the scattering rate 
at T c will result in an increase in spectral weight that 
can more than compensate the decrease due to the open- 
ing of a superconductor gap. We simply adopt the phe- 
nomenological observation of Hosseini et alm& and use the 
following model for the elastic scattering rate, T(T): 



oil:, 



for T >T C 
for T < T r . 



(16) 



Note that we could include a linear temperature depen- 
dence in the normal state, based on the linear resistivity; 
this would serve to increase the (negative) slope in Fig. 2, 
in agreement with observations. However, our attitude 
has been to adopt this simple model to describe the su- 
perconducting state, given that we don't really under- 
stand the normal state, so we will leave this as it is. In 
the superconducting state, one could include a residual 
scattering rate due to impurity scattering, as Hosseini 
et alm& have observed, even in their very clean crystals. 
However, this is a minor effect, and also goes beyond the 
spirit of the simple phenomenology that we are proposing 
here. 

In Fig. 4 we show results for the sum rule with Eq. 
(|16fl "put in by hand". Since we have assumed that the 
crystals are completely free of impurities, the end point 
at zero temperature is the same for all cases shown. It 
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FIG. 3: Plot of sum rule difference, Eq. (|15fl (curves) vs. 
next-nearest neighbour hopping (r = t' /i) for (a) half-filling, 
and (b) n = 0.85. The symbols are for the corresponding dif- 
ference in the magnitude of the kinetic energies. They are in 
good qualitative agreement throughout the range of interest. 



is difficult to try to pin down the amount of inelastic 
scattering that gets suppressed below T c , based on this 
figure, because we can't use the normal state result as a 
benchmark. Our slope is too low in absolute value, and as 
already mentioned, the normal state is undoubtedly not 
reasonably described by our simple approach. Nonethe- 
less, it is clear that the amount of inelastic scattering 
that is suppressed in the superconducting state has to 
exceed a minimum value so that the sum rule will exhibit 
anomalous behaviour. As indicated, a small amount of 



E 



X 
X 

h 

V 



97 - 



95 





i 

S 1/t = 


d-wave 










^s^l/i = 5 meV ^ 




v \\ 
s »\ \ 








— 1/T = 10meV 




i 


^>v^1/T = 15meV 
I 





10000 20000 30000 40000 
T 2 (K 2 ) 
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for various degrees of elastic scattering. Below T c we change 
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temperature is lowered. Further details of the calculation are 
given in the text. 



inelastic scattering to begin with (say, less than 2 meV 
for the parameters used in this figure) , does not make up 
for the decrease expected from conventional BCS theory. 



V. SUMMARY 

It is clear that Fig. 4 "explains" the anomalous optical 
sum rule. To be precise it suggests that the origin of this 
anomaly is the same mechanism that leads to the large 
coherence peak in the microwave^ and to the other ob- 
servations of sharply increased coherence in the supercon- 
ducting stated. It is another matter, however, to claim 
that this increased coherence is an intrinsic part of the 
mechanism. To our knowledge the hole mechanismiiiS 
is the only mechanism for which this is the case in the 
plane. In other models, the increased coherence gener- 
ally arises due to a freezing of the scattering mechanism, 
as originally suggested by Nuss et al. OS Fig. 4 also in- 
dicates a source of non-universal behaviour. If, for some 
materials or doping levels the loss in inelastic scattering is 
small compared to the gap value, then more conventional 
behaviour of the sum rule (eg. the uppermost curve in 
Fig. 4) will prevail. 

The virtue of this work is its simplicity. Without know- 
ing the details of the inelastic scattering process whose 
disappearance in the superconducting state is responsi- 
ble for a sharply decreased scattering rate, we are able 
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to understand the sum rule anomaly. Of course this im- 
mediately implies correlations among the various experi- 
ments. For example, should the anomaly not be present 
in overdoped samples, then the microwave peak and neu- 
tron scattering resonance should also be absent. Im- 
purities should mask the effect, since they presumably 
impact the superconducting state as well as the normal 
stated, so it is important to have very clean samples. We 
have adopted the T 4 dependence for the scattering rate 
suggested in Ref. |23l but the qualitative anomaly will 
not depend on the detailed temperature dependence, but 
rather on the overall magnitude of the change in effective 
scattering rate. 
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